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Although not an intrinsic superconductor, it has been long— known that, when 
intercalated with certain dopants, graphite is capable of exhibiting supercon- 
ductivity [1]. Of the family of graphite— based materials which are known to 
superconduct, perhaps the most well— studied are the alkali metal— graphite in- 
tercalation compounds (GIC) |2] and, of these, the most easily fabricated is 
the CsK system [3] which exhibits a transition temperature T c ~ 0.14 K 
By increasing the alkali metal concentration (through high pressure fabrication 
techniques), the transition temperature has been shown to increase to as much 
as 5 K in C^Na Lately, in an important recent development, Weller et al. 
have shown that, at ambient conditions, the intercalated compounds CqYIo and 
C^Ca exhibit superconductivity with transition temperatures T c ~ 6.5 K and 
11.5 K respectively [5], in excess of that presently reported for other graphite- 
based compounds. We explore the architecture of the states near the Fermi 
level and identify characteristics of the electronic band structure generic to 
GICs. As expected, we find that charge transfer from the intercalant atoms to 
the graphene sheets results in the occupation of the 7r— bands. Yet, remarkably, 
in all those — and only those — compounds that superconduct, we find that an 
interlayer state, which is well separated from the carbon sheets, also becomes 
occupied. We show that the energy of the interlayer band is controlled by a 
combination of its occupancy and the separation between the carbon layers. 

To focus our discussion, we begin with an investigation of the newly discovered CeYb 
system. Our calculations rely on density functional theory (DFT) techniques jij] applied 
within the local density approximation (LDA) Since the 4/-shell of Yb is found to be 
filled (in agreement with experiment), the DFT-LDA can be applied with some confidence. 
Moreover, one may infer that the 4/-orbitals play no essential role in superconductivity. The 
atomic structure of C^Yb, shown in Figured, involves a stacked arrangement of graphene 
sheets intercalated with a triangular lattice configuration of Yb atoms. The stacking of the 
graphene sheets is AAA, while the Yb atoms occupy interlayer sites above the centres of the 
hexagons of the graphene sheets in an ABAB stacking arrangement leading to a P63/MMC 
crystallographic structure. We use the same geometry for C 6 Ca. 

Figure shows the results of the band structure calculation of CgYb centered on the 
Fermi energy. To resolve the qualitative structure of the levels and identify the bands in- 
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FIG. 1: Band structure of CgYb compared to various model systems as inferred from DFT-LDA. 
The panels correspond to (a) atomic structure and (b) geometry of the Brillouin zone of CgYb, band 
structures of (c) expanded Yb, with the lattice parameters and atomic positions corresponding to 
those in CgYb, so that the nearest Yb-Yb distance is increased by 24% as compared to the pure 
FCC structure, (d) the CgYb compound, (e) the "empty" graphite system, where the Yb atoms 
have been removed but the lattice constants are kept fixed, and (f) free electron bands, folded 
into the unit cell of CgYb. The free electron-like band is marked with a thick red line, while pink 
lines indicate hybridized bands with significant charge density in the interlayer region. The band 
structures (except for the free electron bands) are aligned with respect to the core levels and in 
each case the Fermi level is shown with a dashed line. Including correlation effects for the localized 
4/ orbitals within the LDA+U method results only in the expected downward rigid shift of the 
energy of the filled 4/ levels (not shown). 

tersecting the Fermi level, it is helpful to compare the dispersion with the corresponding 
"empty" graphite-like system and the "empty" Yb metal, obtained by simply removing, 
respectively, the Yb and carbon atoms from the CgYb structure whilst holding the lattice 
constants fixed. In CgYb, the previously unoccupied graphite it* bands are seen to inter- 
sect the new Fermi level (due to the enlarged interlayer distance each are almost doubly 
degenerate). Of these, one remains almost unperturbed by the presence of Yb while the 
degeneracy of the other is lifted, and one of the resulting bands strongly hybridizes with a 
new band that also crosses the Fermi level. The latter resembles both the free electron-like 
state of expanded Yb metal (panel(c)) and the interlayer band intrinsic to pure graphite 
(panel (e)). As well as the 5s and 5p orbitals (which lie well outside the range of energies 
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FIG. 2: Band structures of the Li series of GICs, CgCa and MgB2- The latter is shown in a supercell 
that matches the others in the plane of the graphene sheets. The interlayer band is marked with 
a thick red line, while pink lines indicate hybridized bands with significant charge density in the 
interlayer region. The band structures are aligned by the crossing of the tt bands, in each case the 
Fermi level is shown with a dashed line. 

shown in Figure HJi), the 2x74/ orbitals associated with the two Yb atoms in the unit cell 
form a set of localized almost non-dispersive bands which appear at ca. 0.8 eV below the 
Fermi level. Referring to Figure El (panel 4), and the discussion in Ref. jsj], we note that a 
similar phenomenology applies to the band structure of CeCa. 

Although the presence of the interlayer band and its hybridisation with the tc* band helps 
to explain the diminished normal state resistance anisotropy observed in CgYb over that of 
pure graphite , could its occupation be significant for superconductivity? To address this 
question, it is instructive to draw comparison with other superconducting GICs. It is known 



from LEELS measurements of the superconducting compounc 



C«K that, in common with the 



C^Yb system, the interlayer band lies below the Fermi energy jj. However, a more discerning 
test of the significance of the interlayer state is presented by the Li intercalates. While 
intercalation at ambient pressure leads to the CqLi structure, high pressure techniques can 
be used to fabricate CsLi and C2Li compounds. Significantly, while the CsLi system remains 
(like its cousin C@Li) non-superconducting, C2Li exhibits a transition with a T c ~ 1.9 K jioj . 
Referring to Figures|2and|Il one may note that, in common with the CgYb and CgCa system, 
the Li GICs also exhibit a free electron-like interlayer band. The dispersive character of 
this band is reflected clearly in the Fermi surface (see Fig. Ek)- However, significantly, 
increasing Li ion concentration (and, with it, the degree of electron doping) leads to a 
lowering of the interlayer band resulting in its occupancy in the C2Li system. Indeed, 
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FIG. 3: Structure of the interlayer state, (a) Fermi surface of the C2L1 system as inferred from 
the band structure calculation. In its native state, the Fermi surface of a single graphene sheet 
collapses onto two points in the Brillouin zone. Electron-doped, the Fermi surface expands to 
form a set of two cylindrical electron-like surfaces associated with the ir* bands which enclose the 
singular points in the spectrum (red and blue). By contrast, the interlayer band, which adopts a 
nearly free-electron like structure parallel to the planes near the T point, acquires a more complex 
extended character at the Fermi surface (green) . (b) An isosurface of the charge density associated 
with the lowest interlayer band in C2Li at the T point. The colourmap on the left plane shows the 
projection of the band density, with blue corresponding the low and red to high electron density. 
The spatial structure of the interlayer band is essentially identical to that of its unoccupied analogue 
in pure graphite, which is shown in panel (c). 



combining the available data in Table HJ the coincidence of superconductivity in the GICs 
with the occupation of the interlayer state is striking. 

Studies of the free electron-like interlayer state in pure graphite and GICs have a long 
history. Although present in early band structure calculations in the pure system, the 
significance of the (unoccupied) band in graphite and CeLi was first emphasized by Poster- 
nak et al . IllL Il2 | while, subsequently, its existence was confirmed experimentally in both 
graphite [bj and (VLi |bj. (The same band has been equivalently referred to as a metallic 
"sd band" in Ref. [l^J .) Referring to the Li based GICs, it is apparent that the energy, and 



therefore also the occupancy, of this band can vary considerably with respect to the other 
bands present in the above systems. 

It is notable that the dispersive 5d band of Yb in the expanded metal (Fig. ^p) and the 
interlayer state in pure graphite (Fig. |T^) both bear strong resemblance to the dispersion 



5 



GIC 


electron 


interlayer 


interlayer 






doping 


separation 


band occ. 




graphite 





3.35A 


no 
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C 6 Li 
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C 3 Li 
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C 2 Li 


1/2 


3.7A 
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1.9K 
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C 8 K 


1/8 


5.2A 
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0.14K 


C 6 Ca 


1/3 


4.6A 


yes 


11.5K 


C 6 Yb 


1/3 


4.7A 


yes 


6.5K 


C 6 Ba 


1/3 


5.25A 


yes 


? 



TABLE I: Table recording the coincidence of the interlayer state occupation in GICs and the 
phenomenon of superconductivity as inferred from the results of band structure calculation and 
experiment. According to the electronic structure calculations and the observed phenomenology, 
one would predict superconductivity in the GIC CgBa Q]. 

of the band of free electrons in the absence of external potential (Fig. Indeed, we can 
understand the important features of this nearly free-electron like state by studying the 
"empty" graphite system with the intercalant ions removed: In Figure EJ we show the effect 
of charging and layer separation on the energy of the resulting interlayer band. An increase 
in the c axis lattice constant and in the accumulated charge both lead to a lowering of the 
band energy, and thus an increase in its occupancy. It is interesting to note that, for a wide 
range of experimentally known GICs, the interlayer band energy follows closely that of the 
pure graphite model. 

Motivated b y p reliminary electronic structure calculations, previous theoretical studies 
by Al Jishi jl^lisll have emphasized the significance of partially occupied two-dimensional 
graphite it* bands and three-dimensional interlayer bands for superconductivity. Indeed, 
these studies have shown that aspects of superconductivity in some GICs can be described 
well by a two-band phenomenology. Yet the underlying microscopic pairing mechanism 
remains in question. Significantly, unlike the fullerene-based superconducting compounds, 
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FIG. 4: The dependence of the interlayer band energy (at the T point) of the empty graphite 
system with changing c axis spacing for different electron dopings. (Note that the calculation of 
the electron doped system using periodic boundary conditions requires overall charge neutrality, 
enforced by a uniform positive background charge.) Superposed on the figure are data points 
corresponding to different GICs, red and blue crosses indicate whether the compound is found to 
be superconducting or not, respectively. The fractions show the electron doping of the graphite 
sheets in units of electron per carbon atom. In the systems that fall into the shaded area, the 
interlayer band is occupied. 

the 7r* bands are decoupled from the out-of-plane lattice vibrations of the graphene sheets, 
and couple only weakly to the in-plane modes. By contrast, one may expect the interlayer 
band to engage with lattice vibrations of the metal ions and, through their hybridisation with 
the 7T* bands, the graphene sheets. Indeed, recent numerical calculations by Calandra and 
Mauri have revealed the presence of strong electron-phonon coupling in CgCa, facilitated by 
the occupation of the free-electron like band Whether an electron-phonon mechanism 
alone can explain the broad distribution of T c observed across the range of GICs, as well 

n 

as the staging dependence |Tj|, remains a question for future investigation. In this context, 
it is interesting to note that, in those compounds with a low transition temperature, the 
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interlayer band remains nearly orthogonal to the it* bands, as it is in pure graphite, while 
in C^Yb and CeCa the hybridization is strong. 

Intriguingly, through its weak coupling with the graphene layer, one may expect occu- 
pancy of the intrinsic interlayer band of graphite to provide an ideal environment for soft 
charge fluctuations that could promote (s-wave) superconductivity via an excitonic pairing 
mechanism j^. Indeed, the same weak coupling diminishes the potential for charge density 
wave instability which might otherwise compete with the superconducting pair correlations. 
In this context, it is interesting to note that a parallel pairing mechanism involving the ex- 
change of acoustic plasmons was suggested by Takada |21| and was much discussed after the 
discovery of superconductivity in the cuprates. Recently, the model has been re-evaluated for 



layered metals 



22|. Bearing in mind that low- frequency plasma oscillations will, of course, 



couple to the ionic positions — conventional acoustic phonons can be viewed as simply cou- 
pled ionic and electronic plasma oscillations — such processes present an effective mechanism 
to enhance electron-phonon coupling. 

Recognizing the auspicious role played by the free-electron like state in the superconduct- 
ing GICs, it is tempting to look for a common phenomenology in the isoelectronic family of 
boron intercalcates and, in particular, the well-studied superconducting compound MgB 2 

23| . Although the material indeed pos- 



which exhibits a transition temperature T c = 39K 
sesses interlayer states, the band structure (Fig. |2J) shows that it remains unoccupied. How- 
ever, in contrast to CeYb and the other superconducting GICs, in this case charge carriers 
are transferred from the a to the n* band. The resulting band of hole-like a states couples 
strongly to the host lattice providing a natural environment to realise a high temperature 
phonon-mediated superconductivity. 

Finally, it is interesting to note that the interlayer states find their an alog ue as nearly 



free-electron like states in the carbon nanotube system: early theoretical [2J| and experi- 
mental [3] investigations of the alkali-metal nanotube intercalates already indicate a phe- 
nomenology in doping dependence similar to the GICs. On this background, an investigation 
of possible superconductivity in electron-doped carbon nanotubes would seem to be a timely 
and worthwhile enterprise. 

Methods: The results reported here were obtained from an implementation of the LDA 
within the framework of the Castep package 2(|. The methodology, which includes the use 
of ultrasoft pseudopotentials, has been validated in many environments including the rare 
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earth compounds |27(. Here, we take advantage of the fact that, in using ultrasoft pseu- 
dopotentials, one can exercise a considerable degree of freedom in deciding which electron 
shells to treat as fixed core and which to treat as valence, free to participate in bonding. 

The numerical results presented for CgYb were computed treating all the electrons below 
the 5s level for Yb, and below the 2s level for carbon atoms as core, respectively. Plane 
waves up to 480 eV were used to expand the electronic states, and the Brillouin zone was 
sampled with a 12 x 12 x 6 mesh. (A test calculation using a 600 eV cutoff did not yield 



appreciably different results.) To check for possib 
method as implemented in the PWSCF package 



e correlation effects, we used the LDA+U 
28| with U varying between 3 and 7 eV. 



[1] M. S. Dresselhaus and G. Dresselhaus, Intercalation compounds of graphite, Adv. in Phys. 51, 
1-186 (2002). 

[2] Y. Koike, H. Suemetsu, K. Higuchi, and S. Tanuma, Superconductivity in graphite-alkali metal 

intercalation compounds, Physica B+C, 99, 503-508 (1980). 
[3] N. B. Hannay et al, Superconductivity in graphite compounds, Phys. Rev. Lett. 14, 225-226 

(1965). 

[4] I. T. Belash, A. D. Bronnikov, O. V. Zharikov and A. V. Palnichenko, On the superconductivity 
of graphite intercalation compounds with sodium Solid State Commun. 64, 1445-1447 (1987). 

[5] T. Weller, M. Ellerby, S. S. Saxena, R. Smith and N. Skipper, Superconductivity in the inter- 
calated graphite compounds Cq Yband Cq Ca, Nature Physics, to appear. 

[6] M.C. Payne, M. P. Teter, D. C. Allan, T. A. Arias and J. D. Joannopoulos, Iterative mini- 
mization techniques for ab initio total-energy calculations: molecular dynamics and conjugate 
gradients, Rev. Mod. Phys. 64, 1045-1097 (1992). 

[7] J. Perdew and A. Zunger, Self-interaction correction to density-functional approximations for 
many-electron systems, Phys. Rev. B 23, 5048-5079 (1981). 

[8] M. Calandra and F. Mauri, |cond-mat/050 6082 

[9] A. Koma, K. Miki, H. Suematsu, T. Ohno, and H. Kamimura, Density- of- states investigation 
of C8K and occurrence of the interlayer band, Phys. Rev. B 34, 2434-2438 (1986). 
[10] I. T. Belash, A. D. Bronnikov, O. V. Zharikov and A. V. Pal'nichenko, Superconductivity of 
graphite intercalation compound with lithium C^Li, Solid State Commun. 69, 921-923 (1989). 



9 



[11] M. Posternak, A. Baldereschi, A. J. Freeman, E. Wimmer, and M. Weinert, Prediction of 
Electronic Interlayer States in Graphite and Reinterpretation of Alkali Bands in Graphite 
Intercalation Compounds, Phys. Rev. Lett. 50 761-764 (1983). 

[12] N. A. W. Holzwarth, S. G. Louie, and S. Rabii, Interlayer states in graphite and in alkali- 
metal- graphite intercalation compounds, Phys. Rev. B 30, 2219-2222 (1984). 

[13] B. Reihl, J. K. Gimzewski, J. M. Nicholls, and E. Tosatti, Unoccupied electronic states of 
graphite as probed by inverse-photoemission and tunneling spectroscopy, Phys. Rev. B 33, 
5770-5773 (1986). 

[14] Th. Fauster, F. J. Himpsel, J. E. Fischer, and E. W. Plummer, Three- Dimensional Energy 
Band in Graphite and Lithium- Intercalated Graphite, Phys. Rev. Lett. 51, 430-433 (1983). 

[15] S. L. Molodtsov, C. Laubschat, M. Richter, Th. Gantz, and A. M. Shikin, Electronic structure 
of Eu and Yb graphite intercalation compounds, Phys. Rev. B 53, 16621-16630 (1996). 

[16] The structural of C6Ba has been reported in: A. Herold, Crystallo- chemistry of carbon in- 
tercalation compounds, in Intercalated Materials, edited by F. Levy (D. Reidel Publishing 
Company, Dordrecht, Holland, 1979), p. 335; D. Guerard, A. Herold, C.R. Acad. Sc. Paris, 
1974, Serie C, 279, 455. 

[17] R. Al Jishi, Model for superconductivity in graphite intercalation compounds, Phys. Rev. B 
28, 112-116 (1983); 

[18] R. A. Jishi and M. S. Dresselhaus, Superconductivity in graphite intercalation compounds, 

Phys. Rev. B 45, 12465-12469 (1992). 
[19] Y. Iye and S. Tanuma, Superconductivity of graphite intercalation compounds with alkali-metal 

amalgams, Phys. Rev. B 25, 4583-4592 (1982). 
[20] D. Allender, J. Bray and J. Bardeen, Model for an exciton mechanism of superconductivity, 

Phys. Rev. B 7, 1020-1029 (1973). 
[21] Y. Takada, Plasmon Mechanism of Superconductivity in Two- and Three- Dimensional Electron 

Systems, J. Phys. Soc. Jpn. 45, 786-794 (1978). 
[22] For a recent evaluation of the models, see A. Bill, H. Morawitz, and V. Z. Kresin Electronic 

collective modes and superconductivity in layered conductors, Phys. Rev. B 68, 144519 (2003) 
[23] J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani and J. Akimitsu, Superconductivity 

at 39 K in magnesium diboride, Nature (London) 410, 63-64 (2001). 
[24] Y. Miyamoto, A. Rubio, X. Blase, M. L. Cohen and S. G. Louie, Ionic cohesion and electron 

10 



doping of thin carbon tubules with alkali atoms, Phys. Rev. Lett. 74, 2993-2996 (1995). 

[25] H. Shimoda, B. Gao, X. P. Tang, A. Kleinhammes, L. Fleming, Y. Wu and O. Zhou, Lithium 
intercalation into opened single-wall carbon nanotubes: Storage capacity and electronic prop- 
erties, Phys. Rev. Lett. 88, 015502 (2002). 

[26] M. D. Segall, Philip J D Lindan, M J Probert, C J Pickard, P J Hasnip, S J Clark and M C 
Payne, First-principles simulation: ideas, illustrations and the CASTEP code J. Phys. Cond. 
Matt. 14 2717-2744 (2002). 

[27] C. J. Pickard, B. Winkler, R. K. Chen, M. C. Payne, M. H. Lee, J. S. Lin, J. A. White, 
V. Milman and D. Vanderbilt, Structural Properties of Lanthanide and Actinide Compounds 
within the Plane Wave Pseudopotential Approach, Phys. Rev. Lett. 85, 5122-5125 (2000). 



[28] S. Baroni, A. Dal Corso, S. de Gironcoli, and P. Giannozzi, http:/ /www.pwscf.org 



Acknowledgements: We are grateful to S. S. Saxena and Mark Ellerby for bringing 
the CeYb and CeCa compounds to our attention and to Gil Lonzarich for contributing to 
our understanding. 



11 



